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I. INTRODUCTION 

This paper is presented in two parts. The 
first section will be a summary of selected experi- 
mental results obtained on lithium-diffused bulk 
silicon. Particular emphasis was placed on the 
radiation- induced degradation and thermal anneal 
of minority- carrier in bulk silicon because solar 
cell output is related to the minority-carrier life- 
time. The temperature dependence of the 
minority- car rie r lifetime indicates the density and 
energy levels of the recombination centers, pro- 
viding clues to their identity. Electron spin res- 
onance and infrared absorption techniques were 
used to investigate the introduction and anneal of 
three specific radiation- induced defects , which 
are thought to contribute to the recombination 
proces s. 

The results reported herein are a continuation 
and summary of work previously reported (Refs. 
1 - 8 ). 

The second part of the paper will discuss the 
merits of a computer code developed by Gulf 
Radiation Technology which calculates the current- 
voltage (I-V) output of a computer simulated P-N 
junction. With this code, the user stipulates the 
initial properties and environment of a cell. The 
cell's properties and environment may then be 
varied in a manner which might reflect the antic- 


ipated radiation exposure and mission environ- 
ment, and the current- voltage output is recal- 
culated. Thus, a variety of cell types and experi- 
mental conditions can be rapidly evaluated. This 
code can also be utilized by cell manufacturers to 
evaluate the merit of different cell types, doping 
profiles, and cell parameters. It can be utilized 
by those engaged in solar cell testing programs to 
evaluate a greater variety of cell types, radiation 
histories, and ambient conditions than is feasible 
by experimental techniques alone. Systems engi- 
neers can use the code to predict the radiation 
vulnerability of solar cells in varied environments. 

II. EXPERIMENTAL PROGRAM 

A. Samples and Radiation 


The experimental program has been devoted 
to an investigation of bulk lithium- diffused silicon. 
Electron spin resonance and infrared absorption 
were used to study certain specific defects, and 
minority-carrier lifetime measurements were 
used to study radiation- induced recombination 
centers. Samples were lithium- diffus ed by either 
a lithium- oil paint- on technique or by a lithium-tin 
bath technique. Samples with lithium donor den- 
sities from as low as 2 X 10^0 Li/m3 (2 X 1 0 1 4 
Li/cm3) to as high as 4 X 10^3 Li/m^ (4 X 10^ 
Li/cm3) were investigated. Resistivity profiles 
were used to monitor samples for diffusion uni- 
formity. Donor density was occasionally verified 
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by Hall effect measurements. Samples were dif- 
fused from vacuum float- zone (FZ) and quartz- 
crucible (QC) phosphorus-doped silicon, ranging 
in initial room- temperature resistivity from 10 2 
G-m (10^ Sl-cm) to 0. 005 f2-m (0. 5 fl-cm). While 
QC silicon has approximately a hundred times the 
oxygen concentration of FZ silicon, its dislocation 
density is much less. 

Samples were damaged with 30-MeV electrons 
at temperatures ranging from 115 to 300 K or by 
fission neutrons (E > 10 keV) at 273 to 300 K. The 
damage was thermally annealed at temperatures 
up to 673 K. 

B. Study of Three Specific Defects 

The introduction of three specific defects by 
30-MeV electron irradiation and their subsequent 
thermal anneal were investigated using electron 
spin resonance or infrared absorption techniques. 
Specifically, the oxygen- vacancy (Si-Bl), diva- 
cancy (Si-G7), and phosphorus- vacancy (Si-G 8 ) 
centers have been studied. The results of these 
studies have been previously reported elsewhere 
(Refs. 7, 8 ) but may be summarized as follows. 

The oxygen- vacancy and phosphorus- vacancy 
centers were studied by ESR techniques. In QC 
silicon, the Si-Bl center introduction rate is the 
same as in nondiffused silicon (-15 m~ 1 = 0. 15 
cm- 1 ) as long as the lithium density is much less 
than the oxygen density. When the lithium and oxy- 
gen densities are comparable, and B1 center in- 
troduction rate is significantly reduced (- 2 . 5 m" 1 = 
0. 025 cm" 1). This can be attributed to the lithium- 
oxygen pairing, reducing the number of oxygen 
atoms available to form oxygen- vacancy (Si-Bl) 
centers, or to competition between oxygen and the 
positively charged lithium or lithium oxide donors 
for the vacancies. The B1 center in lithium- 
diffused silicon is found to anneal below 400 K in- 
stead of near 600 K, as in nondiffused silicon. 

The LiO + density and conductivity decreased as 
the Si-Bl center annealed. 

Electron spin resonance was also used to study 
the phosphorus- vacancy (Si-G 8 ) center. In 
phosphorus- doped lithium- diffused silicon 
(10 22 P/m 3 = 1 0 1 6 p/ C m3, ~1022 Li/m 3 s 10^6 
Li/cm 3 ), 30-MeV electron irradiations of 10^1 
e/ m^ (1017 e/cm 2 ) below 150 K are found to pro- 
duce Si- G 8 (phosphorus - vacancy ) defects at a rate 
comparable to nondiffused silicon. Production and 
annealing studies on the Si-G 8 center at these den- 
sities and fluence levels indicate that the presence 
of lithium had little or no effect on the creation or 
annealing of that center. It is possible that at 
these high fluence s, relative to the initial lithium 
density, the active lithium has been depleted. 

The divacancy (Si-G7) was observed by moni- 
toring the introduction and anneal of the 1 . 8 - pm 
optical absorption band. In this investigation, it 
was found that for 30-MeV electron irradiations 
at less than 150K to fluences of 10 2 ^ e/m 2 (1017 
e/cm 2 ), the introduction rate of the Si-G7 center 
(divacancy) is comparable in electron- irradiated 
lithium- diffused (5 X 10 22 Li/m 3 = 5 X lO 1 ^ Li/ 
cm 3 ) and nondiffused silicon. The divacancy an- 
neals at or below 300 K in diffused silicon com- 
pared with -325 to 575 K in nondiffused silicon. 

As the 1. 8 - pm divacancy band disappears, new 


bands near 1. 4 and 1. 65 pm appear, and these 
anneal near 600 K. These bands have been ob- 
served by other investigators (Ref. 9) and attri- 
buted to the divacancy plus one or two lithium 
atoms. 

C. Minority- Car rier Lifetime Measurements 

The introduction and anneal of recombination 
centers were determined from minority- carrier 
lifetime measurements. Minority- carrier lifetime 
was determined from measurements of the decay 
of photoconductivity or by steady- state techniques. 
These measurements are described in detail else- 
where (Refs. 3, 8 ). The present results were at 
injection levels (An/ng) of less than 1%. Inter- 
pretation of the results is based on Shockley- Read 
recombination theory (Ref. 10). Capture cross- 
section temperature dependence as estimated by 
Lax (Ref. 11) or Leadon (Ref. 12) was assumed. 
The temperature dependence of minority-carrier 
lifetime was measured before and after irradiation 
and again after thermal anneal. Figure 1 shows 
minority-carrier lifetime versus inverse tem- 
perature for a lithium- diffused quartz- c rucible 
and a lithium- diffused float- zone sample. The 
minority-carrier lifetime temperature dependence 
of the QC silicon can be attributed to recombina- 
tion centers further than about 0. 3 eV from either 
band edge, but the lifetime of FZ silicon has a 
more severe temperature dependence and a dis- 
continuity near 1000/T = 4.5K’ 1 . This can be 
attributed to recombination centers about 0. 17 eV 
below the conduction band, which were not ob- 
served in QC silicon, in addition to centers fur- 
ther than 0. 35 eV from either band edge. The 
lines drawn through the data represent the low 
injection-level lifetimes calculated by Shockley- 
Read theory, assuming the temperature- dependent 
capture cross sections of Lax. The presence of 
a recombination center near E c - 0. 17 eV in FZ 
silicon and its absence in QC silicon suggests that 
this center is not the oxygen- vacancy (Si-Bl) cen- 
ter, which is known to have an energy level in this 
range. 

The minority-carrier lifetime degradation 
rate is defined in terms of lifetime T and electron 
fluence tj> to be 


compared with phosphorus-doped silicon, and was 
found to increase with increasing lithium donor 
density. The degradation constant was found to be 
essentially independent of oxygen content. This 
shows that the 30-MeV electron- induced defects 
either contain lithium or are affected in their pro- 
duction by lithium. The minority- car rier lifetime 
degradation constant was found to vary with tem- 
perature in a way which could be attributed to the 
temperature dependence of the capture cross 
section. 


Both isothermal and isochronal annealing ex- 
periments have been performed on a variety of 
lithium- diffused samples. The unannealed fraction 
of annealable defects for isothermal annealing is 
given by 
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where tq is the lifetime after irradiation and be- 
fore anneal, tj is the fully annealed lifetime, and 
T t is the lifetime at the annealing temperature after 
time t. Linearity of the isothermal annealing data 
implies first-order annealing kinetics (Ref. 13) 
in which the number of recombination centers is 
given by 


A more detailed discussion of our experimen- 
tal findings and theoretical interpretations is pre- 
sented elsewhere (Ref. 8), and a considerable body 
of additional data, provided by us and other inves- 
tigators using other techniques, is available on 
bulk silicon. The following section of this paper 
describes one means of applying this body of data 
to the problem of predicting the response of a 
solar cell to various radiation and environmental 
conditions . 

III. CODE FOR PREDICTING PERFORMANCE 
OF SOLAR CELLS 


N = N Q exp (-Rt) 


where Nq is the initial number of annealable defects 
and R is the rate constant, which is given by 


R = v exp (-E/kT) 


where E is the activation energy and v is the effec- 
tive frequency factor. In every case, our iso- 
thermal anneals indicated first-order annealing 
kinetics were appropriate. 

For isochronal anneals, the unannealed frac- 
tion of annealable defects is given by 


1/t 

177 


T. 

0 


1/t 


1/t 


where is the lifetime observed at 300 K after 
anneal at temperature T, and tq and are the 
pre- and post-anneal lifetimes. 

Analysis of isochronal and isothermal anneals, 
using first-order kinetics, yields the activation 
energy and frequency factor of the rate constant. 
Figure 2 shows the results of five separate iso- 
chronal annealing experiments on lithium- diffused 
QC silicon. A single first-order annealing stage 
at 380 ±10 K is readily apparent. Similar results 
were obtained on fission- neutron- irradiated 
silicon. 

Table 1 summarizes our experimental findings 
for the annealing of radiation- induced recombina- 
tion centers in lithium- diffused silicon exposed to 
30-MeV electron and fission neutrons (E > 10 keV). 
The activation energies compare with the 0. 65 
± 0. 05 eV energy for lithium diffusion in oxygen- 
lean FZ silicon, and with 1. 07 ± 0. 05 eV energy of 
lithium- oxygen dissociation and diffusion in oxygen- 
rich QC silicon (Refs. 14-16). This data supports 
the supposition that radiation- induced damage, in- 
cluding the clustered damage attributed to fission 
neutrons, is annealed by the diffusion of lithium to 
the defect. 

In addition, the effective frequency factors are 
much less than the atomic frequency factor of about 
10 13 s - 1 and scale with the lithium donor density. 
This is consistent with a process involving the 
long-range migration of one annealing species to 
another (Ref. 13), which we take to be the migra- 
tion of lithium to the damage site. 


The purpose of this presentation is to de- 
scribe a computer program that can be used for 
predicting the output performance of solar cells. 
This code was developed under Defense Atomic 
Support Agency contract and applied to solar cells 
under the JPL contract. Briefly, this code, which 
is fully operational at the present time, solves the 
complete equations for a one-dimensional P-N 
junction. The code starts with the current and 
continuity equations for electrons and holes, 
Poisson's equation for the electric field inside the 
cell, and appropriate equations for the generation 
and recombination of excess carriers. Arbitrary 
one-dimensional space profiles of doping density, 
mobilities, recombination lifetimes, and ioniza- 
tion densities are inputs to the code. The partial 
differential equations are converted to finite dif- 
ference approximations, which are then solved by 
iteration. In principle, there are no adjustable 
parameters in the code. If one had accurate values 
for all of the material properties, the code would 
predict the absolute magnitudes of the current and 
voltage for selected values of load resistance. 

Thus, the accuracy of the results is primarily 
determined by how well one can measure or esti- 
mate the important physical parameters of the 
cell. 


This code could be useful both in the design of 
new devices and in the analysis of radiation effects 
to solar cells. For example, if a cell manufac- 
turer is designing a new type of cell, perhaps with 
an unusual doping profile, he could estimate the 
performance of the new device relative to a con- 
ventional cell and optimize its parameters before 
he goes through the expense of constructing it. 
Similarly, the radiation- effects physicist could 
use such a code to help him analyze his solar cell 
data in terms of kinds and densities of defects 
created by the radiation and to extrapolate the re- 
sults to other radiation environments. 

The three basic equations of the code are the 
one- dimensional continuity equations for holes (p) 
in the valence band and electrons (n) in the con- 
duction band, and Poisson's equation for the elec- 
tric field (E). 


9p(x, t) 

g (x, t) - R(x, t) - 

9J (x, t) 
P 

(i) 

at 

3x 

9n(x, t) 

g(x, t) - R(x, t) - 

3J (x, t) 
n 

(2) 

at 

3x 

9E(x, t) 
dx 

|L[p(x, t) - n(x, t) 

+ AN(x)] 

(3) 
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In the above equations, the dimension x is normal 
to the face of the cell and g(x, t) is the time- and 
space- dependent ionization rate of electron-hole 
pairs by the external radiation. This quantity is 
calculated ahead of time, using the spectral dis- 
tribution of the incident light and the attenuation 
coefficients for the particular material, and is 
then an input into the code. The quantity q is the 
magnitude of the electronic charge, K is the di- 
electric constant in MKS units, AN(x) is the net 
doping of the semiconductor (positive for n-type, 
negative for p-type), and Jp and J n are the particle 
current densities given by 


J 


P 


pp. E - 
P 


D 35 
p3x 


(4) 


J 

n 


-nu E 
"n 


- D I 2 

n ox 


(5) 


In Eqs. (4) and (5), the first terms are the drift 
components of the currents with mobilities (ip and 
p n , while the second terms give the diffusion cur- 
rents with diffusion coefficients Dp and D^. 

The term R(x, t) is the rate of recombination 
of excess carriers via recombination centers. 
There are several different options available to a 
user, depending on the sophistication he wishes to 
use in describing the recombination process. In 
many situations, it is adequate to use a Shockley- 
Read-type recombination rate for a single recom- 
bination level:. 


2 

np - n. 

R " t q (P + p F ) + t q (n + n ) (6) 

n p 

where n^ is the intrinsic carrier density, pjr and 
np are the values of p and n when the Fermi level 
coincides with the level of the recombination cen- 
ter that is being considered, and tq and tq are 
the recombination lifetimes of electrons in hlavily 
p-type material and of holes in heavily n-type 
material, respectively. 

This recombination equation does have the 
limitation that it only simulates the simultaneous 
annihilation of a free electron and a free hole. It 
does not consider the trapping of carriers on de- 
fect centers and, therefore, carrier removal. 
Since the latter phenomena can be important when 
the details of the trapping kinetics affect the car- 
rier lifetimes or removal rates, the code also has 
the capability of simulating both single-level and 
double- level traps. This could be important, for 
example, if the radiation produces a high-density 
region of defect centers which significantly de- 
pletes the majority carriers and creates a p-n 
junction effect. 

For a single- level trap with density Nj and 
two possible charge states, negative (Nj) and 
neutral (N°), the (-R) term in Eq. (1) would be 
replaced by 


- R P = %( pN i ' n i p f) (7) 


and the (-R) term in Eq. (2) would be replaced by 


-R 


o o / 

= a i 
n n\ 


nlNfj - N 


l n F) 


(8) 


The quantity a' is the product of the thermal ve- 
locity of a hole in the valance band and the cross 
section o' for capture of a hole by a negatively 
charged defect, and o° is a similar term for cap- 
ture of an electron by a neutral defect. In addition 
to modifying Eqs. (1) and (2), the density of the 
negatively charged center (N j ) must be included 
in the summation of charges in Eq. (3). Further- 
more, additional equations are required for the 
rate of change of Nj and N?. For constant total 
Nj, 


dN , dN° 

L _ L _ R ° 

dt dt _ n 


(9) 


Similar equations are available in the code for a 
double-level defect with three charge states, N 2 , 
N°, and Nj, but, for brevity, they are not pre- 
sented here. 


To obtain a solution to the above equations, 
they are converted to finite difference form for 
some prescribed mesh distribution. The distri- 
bution of the meshes can be chosen arbitrarily by 
the user to most judiciously utilize the permissible 
number of meshes. On the Univac 1108, core 
storage limits the number of meshes to about 230, 
which is an adequate number to simulate a three- 
layer planar transistor, and is certainly enough 
for a two-region solar cell. 

Since this code was originally developed to 
study time- dependent processes in semiconductor 
devices, the finite difference equations are solved 
by iteration for finite time steps. Under condi- 
tions of steady illumination and circuit load, the 
code quickly goes to the steady- state solution. 

The resulting steady voltage and current form one 
point on the I- V curve. Additional points are ob- 
tained by changing the external load resistance 
and again allowing the problem to run to equilib- 
rium. 

Since Eqs. (1), (2), and (3) are time- 
dependent partial differential equations, one must 
specify boundary conditions at the two ends of the 
device and self-consistent initial distributions on 
n, p, and E. The self-consistent initial distri- 
butions could have been specified in any number of 
ways. To simplify the task of the user, a new 
problem is started from a condition of zero elec- 
tric field and charge neutrality everywhere. This 
is usually a physically unrealistic, but mathe- 
matically correct, starting point. The code then 
does the work of proceeding to the correct solu- 
tion. For a subsequent run, with the same device, 
it is usually possible and economical to restart 
from any previous cycle. 

Two types of boundary conditions are avail- 
able. The first, called "bulk, " forces the car- 
rier densities at the two ends of the device to have 
zero slopes. This is the condition that would 
apply if the p-n junction were far from the contacts 
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and the densities have a chance to approach their 
bulk values far from the junction (hence the name 
"bulk"). This boundary condition is suitable when 
the details of the contact are not important in the 
problem. 

On the other hand, when the details of the 
contact between the semiconductor and the metal 
are important, the contact potential between the 
two materials can be simulated in the second type 
of boundary condition. A large surface recom- 
bination velocity can be simulated by defining a 
narrow region at the surface with a very short re- 
combination lifetime. Of course, an additional 
requirement is that the current through the ex- 
terior circuit must equal the current through the 
device, including the displacement current in time- 
dependent problems. Detailed descriptions of the 
starting and operating procedures for the code can 
be found in Refs. 17 and 18. 

Since this code solves the complete equations 
for the solar cell, the answers that it provides 
should be as good, within the accuracy of the finite 
difference method of solution, as the input data 
supplied to the code. For comparative design 
studies, the theoretical device characteristics can 
be specified exactly and the effect of varying one 
or more parameters individually can be evaluated. 
For the analysis of radiation effects in an actual 
device, the problem is somewhat different because 
the parameters for the undamaged test device, 
such as doping profile and lifetimes, are usually 
not accurately known. For the best results, it is 
always desirable to avail oneself of all possible 
nondestructive tests to determine the device 
parameters, such as capacitance- voltage mea- 
surements and resistivity for the doping density 
profile and diffusion lengths for the lifetimes. 

However, in the final analysis, it may be 
necessary to adjust these parameters somewhat 
by trial and error to obtain acceptable agreement 
with the preirradiated test data. Once these char- 
acteristics are determined for the undamaged cell, 
the change in the bulk parameters due to radiation 
damage can then be estimated and included in the 
code. If the test data on radiation effects in the 
bulk material are sufficiently good so that the 
model for the radiation damage and the change per 
unit dose can be correctly determined, the pre- 
dicted results for the irradiated cell should be in 
acceptable agreement with the post-irradiation test 
data. If there is significant disagreement, it must 
mean that the model for the radiation damage or 
the magnitude of the radiation effects is different 
from that estimated. In this way, the code can be 
used as a tool for analyzing and interpreting radia- 
tion effects test data. After agreement has been 
obtained between the experimental and calculated 
results for a number of test conditions, the code 
can then be used to extend the test results to other 
test or operating conditions. 

As an example of the results from this code. 
Fig. 3 shows predicted I-V curves for a particular 
solar cell for several fluences of 300-keV protons 
impinging on the illuminated surface. The param- 
eters for the undamaged cell were adjusted to give 
agreement with the performance curve supplied by 
the manufacturer. The proton damage was simu- 
lated by a distribution of recombination centers, 
concentrated mainly in the region of the maximum 


range of the protons. If adequate test data for the 
proton fluences were available, the validity of the 
model for the proton damage could be tested by 
comparing experiment to prediction. 
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Table 1. Summary of the results of isochronal and isothermal annealing experiments on electron- and 
neutron-irradiated lithium-diffused float-zone and quartz -crucible silicon 


Sample type 

Float zone 

Quartz crucible 

Lithium donor 

2. 5 X 10 22 

4. 5 X 10 20 

2 X 10 21 

2.5 X 10 21 

4. 5 X 10 20 

2. 5 X 10 21 

density, Li/m 







Fluence 
(particles/m ) 

~io 17 

10 16 to 10 17 

5 X 10 14 

1. 5 X 10 16 
to 8 X 10 16 

io 16 

14 

4 X 10 

Radiation 

30-MeV 

30-MeV 

Fission 

30-MeV. 

30-MeV 

Fission 


electrons 

electr ons 

neutrons 

electrons 

electrons 

neutrons 

Activation 

0. 85 ±0. 10 

0. 60 ±0. 10 

0. 67 ±0. 05 

1. 0 ±0. 2 

0. 75 ±0. 10 

1. 2 ±0. 6 

energy, eV 







Frequency factor, 
-1 

sec 

£ 

o 

o 

Unknown 

~10 7 

10 8 to 10 10 

10 7 to 10 8 

Unknown 
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MINORITY-CARRIER LIFETIME 





UNANNEALED FRACTION OF ANNEALABLE DEFECTS 



360 380 

ANNEAL TEMPERATURE (°K) 


Unannealed fraction of annealable defects after 5-min isochronal anneals at indicated 
temperatures after 30-MeV irradiation of lithium- diffused N-type QC silicon 
irradiated at 300 K 





CURRENT (A/cm 2 ) 


Fig. 
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